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Abstract. This chapter is concerned with the question: how do gravitational waves 
(GWs) interact with their detectors? It is intended to be a theoretical review of the 
fundamental concepts involved in interferometric and acoustic (Weber bar) GW an- 
tennas. In particular, the type of signal the GW deposits in the detector in each case 
will be assessed, as well as its intensity and deconvolution. Brief reference will also be 
made to detector sensitivity characterisation, including very summary data on current 
state of the art GW detectors. 



1 Introduction 

Gravitational waves (GW), on very general grounds, appear to be a largely un- 
avoidable consequence of a well established fact that no known interaction prop- 
agates instantly from source to observer: gravitation would be the first exception 
to this rule, should it be described by Newton's theory. Indeed, the Newtonian 
gravitational potential <f>(x,t) satisfies Poisson's equation 

V 2 (f>(x,t) = -4nGg(x,t) , (I) 

where g(x,t) is the density of matter in the sources of gravitational fields, and 
G is Newton's constant. But, since (^) contains no time derivatives, the time 
dependence of cj>(x,t) is purely parametric, i.e., time variations in g(x,t) in- 
stantly carry over to cf>(x,t), irrespective of the value of x. So, for example, 
non-spherically symmetric fluctuations in the mass distribution of the Sun (such 
as e.g. those caused by solar storms) would instantly and simultaneously be felt 
both in the nearby Mercury and in the remote Pluto. . . 

Quite independently of the quantitative relevance of such instant propagation 
effect in this particular example - which is none in practice -, its very existence 
is conceptually distressing. In addition, the asymmetry between the space and 
time variables in (|l|) does not even comply with the basic requirements of Special 
Relativity. 

Einstein's solution to the problem of gravity, General Relativity (GR), does 
indeed predict the existence of radiation of gravitational waves. As early as 1918, 
Einstein himself provided a full description of the polarisation and propagation 
properties of weak GWs [[[] . According to GR, GWs travel across otherwise flat 
empty space at the speed of light, and have two independent and transverse 
polarisation amplitudes, often denoted h + {x,t) and h x (x,t), respectively 
In a more general framework of so called metric theories of gravity, GWs are 
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allowed to have up to a maximum 6 amplitudes, some of them transverse and 
some longitudinal M . 

The theoretically predicted existence of GWs poses of course the experimen- 
tal challenge to measure them. Historically, it took a long while even to attempt 
the construction of a gravitational telescope: it was not until the decade of the 
1960's that J. Weber first took up the initiative, and developed the first grav- 
itational antennas. These were elastic cylinders of aluminum, most sensitive to 
short bursts (a few milliseconds) of GWs. After analysing the data generated 
by two independent instruments, and looking for events in coincidence in both, 
he reported evidence that a considerably large number of GW flares had been 
sighted [|. 

Even though Weber never gave up his claims of real GW detection ||, his 
contentions eventually proved untenable. For example, the rate and intensity of 
the reported events would imply the happening of several supernova explosions 
per week in our galaxy |(| , which is astrophysically very unlikely. 

It became clear that more sensitive detectors were necessary, whose design 
and development began shortly afterwards. In the mid 1970's and early eighties, 
the new concept interferometric GW detector started to develop which 
would later lead to the larger LIGO and VIRGO projects, as well as others 
of more reduced dimensions (GEO-600 and TAMA), and to the future space 
antenna LISA. In parallel, cryogenic resonant detectors were designed and con- 
structed in several laboratories, and towards mid 1990's the next generation of 
ultracryogenic antennas, NAUTILUS (Rome), AURIGA (Padua), ALLEGRO 
(Baton Rouge, Louisiana) and NIOBE (Australia), began taking data. More re- 
cently, data exchange protocols have been signed up for multiple detector coinci- 
dence analysis ||. Based on analogous physical principles, new generation spher- 
ical GW detectors are being programmed in Brazil, Holland and Italy (|lO|Jll|Jl^] . 

In spite of many years of endeavours and hard work, GWs have proved elusive 
to all dedicated detectors constructed so far. However, the discovery of the binary 
pulsar PSR 1913+16 by R. Hulse and J. Taylor in 1974 @, and the subsequent 
long term detailed monitoring of its orbital motion, brought a breeze of fresh air 
into GW science: the measured decay of the orbital period of the binary system 
due to gravitational bremsstrahlung accurately conforms to the predictions of 
General Relativity. Hulse and Taylor were awarded the 1993 Nobel Prize in 
Physics for their remarkable work. As of 1994 [jl4j, the accumulated binary pulsar 
data confirm GR to a high precision of a tenth of a percent]]. 

The binary pulsar certainly provides the most compelling evidence to date 
of the GW phenomenon as such, yet it does so thanks to the observation of 
a back action effect on the source. Even though I do not consider accurate the 
statement, at times made by various people, that this is only indirect evidence of 
GWs, it is definitely a matter of fact that there is more to GWs than revealed by 
the binary pulsar. . . For example, amplitude, phase and polarisation parameters 

1 It appears that priorities in pulsar observations have since shifted to other top- 
ics of astrophysical interest, so it is difficult to find more recent information on 
PSR 1913+16. 
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of a GW can only be measured, according to current lore, with dedicated GW 
antennas. 

But how do GW telescopes interact with the radiation they are supposed to 
detect? This is of course a fundamental question, and is also the subject of the 
present contribution, where I intend to review the theoretical foundations of this 
problem, and its solutions as presently understood. In section ^, I summarise the 
essential properties of GWs within a rather large class of possible theories of the 
gravitational interaction; section || briefly bridges the way to sections | and | 
where interferometric and acoustic detector concepts are respectively analysed 
in some detail. Section || also includes aboundant reference to new generation 
spherical detectors in its various variants (solid, hollow, dual). For the sake of 
completeness, I have added a section (section ||) with a very short summary of 
detector characterisation concepts, so that the interested reader gets a flavour 
of how sensitivities are defined, what do they express and how do GW signals 
compare with local noise disturbances in currently conceived detectors. Section [7] 
closes the paper with a few general remarks. 



2 The nature of gravitational waves 

Quite generally, a time varying mass-energy distribution creates in its surround- 
ings a time varying gravitational field (curvature). As already stressed in sec- 
tion [l], we do not expect such time variations to travel instantly to distant places, 
but rather that they travel as gravity waves across the intervening space. 

Now, how are these waves "seen"? A single observer may of course not feel 
any variations in the gravitational field where he/she is immersed, if he/she is 
in free fall in that field - this is a consequence of the Equivalence Principle ]l5| . 
Two nearby observers have instead the capacity to do so: for, both being in free 
fall, they can take each other as a reference to measure any relative accelerations 
caused by a non-uniform gravitational field, in particular those caused by a 
gravitational wave field. We can rephrase this argument saying that gravitational 
waves show up as local tides, or gradients of the local gravitational field at the 
observatory. 

In the language of Differential Geometry, tides are identified as geodesic de- 
viations, i.e., variations in the four- vector connecting nearby geodesic lines. It is 
shown in textbooks, e.g. 0], that the geodesic deviation equation is 

^2 — ^iivpa X X ^ , (2) 

where "Z?" means covariant derivative, s is proper time for either geodesic, -R M „ p(T 
is the Riemann tensor, x v is a unit tangent vector (again to either geodesic), 
and £ M is the vector connecting corresponding points of the two geodesic lines. 

The GW fields considered in this paper will be restricted to a class of per- 
turbations of the geometry of otherwise flat space-time, with the additional as- 
sumptions that they be 
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- small 

- time-dependent 

- vacuum perturbations 



This is certainly not the most general definition of a GW yet it will suffice 
to our purposes here: any GWs generated in astrophysical sources and reaching 
a man-made detector definitely satisfy the above requirements. The interested 
reader is referred to Jl(| for a thorough treatment of more general GWs. 

Following the above assumptions, a GW can be described by perturbations 
h^i,(x,t) of a flat Lorentzian metric rj^, i.e., there exist quasi- Lorentzian coor- 
dinates (x,t) in which the space-time metric g^^x^t) can be written 

g, iV {x,t) = r} IJ , v + h llv (x,t) (3) 

with 



rip* = diag (-1,1,1,1) , \h^{x,t)\ < 1 . (4) 

The actual effect of a GW on a pair of test particles is, according to (||), 
determined by the Riemann tensor R pvpa , and this in turn is determined by the 
functions h pv . I now review briefly the different possibilities in terms of which 
is the theory underlying the physics of gravity waves, i.e., which are the field 
equations which the h pv satisfy. 

2.1 Plane GWs according to General Relativity 

The vacuum field equations of General Relativity are, as is well known 

R»» = , (5) 

where R pl/ is the Ricci tensor of the metric g^. If quadratic and successively 
higher order terms in the perturbations are neglected then this tensor can 
be seen to be given by 

= Dh^ - dpdji^ - d^djfp , (6) 

with 

□ = rf d„d v , V = V ~ihri^, h = r,^ . (7) 
New coordinates (x' , t') can be defined by means of transformation equations 



x'* 1 =x"+e^{x,t) , (8) 

and these will still be quasi-Lorentzian if the functions e^{x,t) are sufficiently 
small. More precisely, the GW components are, in the new coordinates, 
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provided higher order terms in e M are neglected. Thus "sufficiently small" means 
that the derivatives of be of the order of magnitude of the metric perturbations 



hfj, u ~ so that in the new coordinates ; 
9uu = Vi^v + h'^. It is now possible, see 
way that the gauge conditions 



the metric tensor also splits up as 
to choose new coordinates in such a 



cUp = 



hold. This being the case, equations (||) read 







(10) 



(11) 



which are vacuum wave equations. Therefore GWs travel across empty space 
at the speed of light, according to GR theory. Plane wave solutions to ( [TT| ) 
satisfying ( |Io| ) can now be constructed M which take the form 



hfj(x,t) 



/0 







h+(t- 
h x (t- 



z) 
z) 





-h+(t 
o 



z) 
-z) 







(12) 





\0 0/ 

for waves travelling down the z-axis. The label TT is an acronym for transverse- 
traceless, the usual denomination for this particular gauge. 

The physical meaning of the polarisation amplitudes in ( |l2] ) is clarified by 
looking at the effect of an incoming wave on test particles. Consider e.g. two 
equal test masses whose center of mass is at the origin of TT coordinates; let 
4) be their distance in the absence of GWs, and (9, <p) the orientation (relative 
to the TT axes) of the vector joining both masses. Making use of the geodesic 
deviation equation (0), with the Riemann tensor associated to (|l~2"|), it can be 
seen that the GW only affects the transverse projection of the distance relative 
to the wave propagation direction (the z-axis); in fact, if £{t) = (£ M S,^) 1 ^ 2 is such 
distance then some simple algebra leads to the resultn 



t{t) = l 



1 



{h + (t) cos 2ip + h x (t) sin 2^} sin 2 9 



(13) 



It is very important at this point to stress that the wavelength A of the 
incoming GW must be much larger than the distance £(t) between the particles 
for (O) to hold, i.e., 



t{t) < A , 



(14) 



and this is a condition which must be added to the already made assumption 
that |/i M „| < 1. 

A graphical representation of the result ( |l3| ) is displayed in figure [j]: a num- 
ber of test particles are evenly distributed around a circle perpendicular to the 

2 Note that the Riemann tensor is calculated at the center of mass of the test particles, 
therefore at x = 0. But it can also be calculated at the position of either one of 
them - this would only make up for a negligible second order difference. 
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Fig. 1. The "+" and "x" polarisation modes of a GW, according to GR theory. The 
dotted lines (circles) indicate the test particles' position in the absence of GW signal. 
Each step in the graph corresponds to a quarter of the period of the driving GW, as 
labeled atop. 

incoming GW, i.e., in the xy plane. When a periodic signal comes in, the dis- 
tances between those particles change following (|l3|); note that the changes are 
modulated by the angular factors, i.e., according to the particles' positions on 
the circle. The "+" mode is characterised by a vanishing wave amplitude h x , 
while in the "x" mode h+ vanishes. 

2.2 Plane GWs according to metric theories of gravity 

Although General Relativity has never been questioned so far by experiment, 
there are in fact alternative theories, e.g. Brans-Dicke theory |l7), which are 
interesting for a number of reasons, for example cosmological reasons jl8| . Gen- 
erally, these theories make their own specific predictions about GWs, and they 
partly differ from those of GR - just discussed. The term metric theory indi- 
cates that the gravitational interaction affects the geometry of space-time, i.e., 
the metric tensor g^ is a fundamental ingredient - though other fields may also 
be necessary to complete the theoretical scheme. Obviously, General Relativity 
falls within this class of theories. 

The appropriate scheme to assess the physics of such more general class of 
GWs was provided long ago jL9| . The idea is to consider only plane gravitational 
waves, which should be an extremely good approximation for astrophysics, given 
our great distance even to the nearest conceivable GW source, and to characterise 
them by their Newman- Penrose scalars p(| . 

It appears that only six components of the Riemann tensor out of the usual 20 
are independent in a plane GW; these are given by the four Newman-Penrose 
scalars 
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V2(v),y 3 (v),y 4 (v) y <t> 22 (v) , (15) 

of which y V 2 and 4> 2 2 are real, while SP3 and ¥4 are complex functions of the 
null variable v - sec ^(J for all notation details. If a quasi-Lorentzian coordinate 
system is chosen such that GWs travel along the z-axis then v = t — z, and one 
can calculate the scalars (|15l) to obtain 



¥ 2 (t- 


z) = 


g Rfztz 1 


*s(t- 


*) = 


~\ (~Rtxtz + ' l Rtytz) •> 




z) = 


— Rtxtx + ^iyiy + ^-Rtxty 




*) = 


Rtxtx Rtyty • 



(16a) 
(16b) 
(16c) 
(16d) 

General Relativity is characterised by — z) being the only non- vanishing 
scalar, while in Brans-Dicke theory 4> 2 2(£ — z) also is different from zero - see 
for full details. 

It is relevant to remark at this stage that the only non-trivial components 
of the Riemann tensor of a plane GW are the so called "eZeciric" components, 



Rtitj, as we see in equations (16a-d) above. These are, incidentally, also the 
only ones which appear in the geodesic deviation equation (|2|), since one may 
naturally choose = (1; 0, 0, 0) (]. 

This fact helps us make a graphical representation of all six possible polar- 
isation states of a general metric GW in the same manner as in figure [l]. The 
result is represented in figure || - whose source is reference [|| . The idea is to 
take a ring of test particles, let a GW pass by, and analyse the results of the 
displacements it causes in the distributions of those particles, just as done in 



section 2.2. Note that the first three modes are transverse, while the other three 
are longitudinal - see the caption to the figure. As already stressed, GR only 
gives rise to the two ¥4 modes. 



3 Gravitational wave detection concepts 

We are now ready to discuss the objectives and procedures to detect GWs: 
knowing their physical structure, one can design systems whose interaction with 
the GWs be sufficiently well understood and under control; suitable monitoring 
of the dynamics of such systems will be the source of information on whatever 
GW parameters may show up in a given observation experiment. 

There are two major detection concepts: interferometric and acoustic detec- 
tion. Historically, the latter came first through the pioneering work of J. Weber, 
but interferometric GW antennas are at present attracting the larger stake of 

J Note that, with this choice, R^ pv(T x^x" = Rt P ta but, because of the symmetries of 
the Riemann tensor, only values of p and a different from zero, i.e., {per} = {ij}, 
give non-zero contributions. 



8 J. Alberto Lobo 







Fig. 2. The six polarisation modes of a plane, metric GW. The dotted and solid lines 
correspond to distributions separated by a half period of the incoming wave, and it is 
assumed that the particles lie on a circle when no GW is acting. Note the indicators 
that the modes in the upper row are excited by a wave which arrives perpendicular 
to the particles' plane, while those in the lower row correspond to GWs which are 
in the same plane as the particles. It must however be clarified that the mode <J>22 
has spherical symmetry, so it includes a combination of longitudinal and transverse 
excitations in like proportions. 



the investment in this research field, both in hardware and in human commit- 
ment. This is because much hope has been deposited in their capabilities to reach 
sufficient sensitivity to measure GWs for the first time. 

Interferometric detectors aim to measure phase shifts between light rays 
shone along two different (straight) lines, whose ends are defined by freely sus- 
pended test masses. This is done in a Michelson layout, using mirrors, beam 
splitters and photodiodes. Acoustic detectors are instead based on elastically 
linked test masses - rather than freely suspended - which resonantly respond to 
GW excitations. 

These qualitative ideas can be made quantitatively precise, but the process 
is a non-trivial one and has important subtleties which must be properly un- 
derstood for a thorough assessment of the detector workings and readout. The 
next sections are devoted to explain with some detail which are the theoretical 
principles governing the behaviour and response of both kinds of GW antennas. 
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Fig. 3. The "naive" interferometric detector concept: a GW coming perpendicular to 
the sheet's plane, "+" polarised relative to the x and y axes, causes the end masses Mi 
and Mi to oscillate in phase opposition relative to the central mass Mo. Light is shone 
into the system, and suitable beam splitters and mirrors are attached to the masses; 
length changes are then measured interferometrically, which directly lead to determine 
the GW amplitude. 

4 Interferometric GW detectors 

A rather naive idea to measure the effect of an incoming GW is provided by 
the following argument - see also figure |^ for reference: let a GW having a "+" 
polarisation (assume GR for simplicity at this stage) come in perpendicular to 
the local horizontal at a given observatory; if three masses are laid down on the 
vertices of an ideally oriented isosceles right triangle then, as we saw (figure |l|), 
the catets shrink and stretch with opposite phases. If a beam of laser light is 
now shone into the system, and a beam-splitter attached to mass Mq and mirrors 
attached to Mi and M2, then one can think of measuring the distance changes 
between the masses by simple interferometry. 

This may look like a very reasonable proposal for a detector yet the follow- 
ing criticism readily suggests itself: gravitation is concerned with geometry, i.e., 
gravitational fields alter lengths and angles; therefore GWs will affect identi- 
cally both distances between the masses of figure ^ and the wavelength of the 
light travelling between them thus leading to a cancellation of the conjectured 
interferometric effect. . . 

While the criticism is certainly correct, the conclusion is not. The reason is 
that it overlooks the fact that gravitation is concerned with the geometry of 
space-time - not just space. In the case of the above GW it so happens that, 
in TT coordinates, the time dimension of space-time is not warped in the GW 
geometry - see the form of the metric tensor in equation ( |l^ ) - while the trans- 
verse space dimensions are. Consequently an electromagnetic wave travelling in 
the xy plane experiences wavelength changes depending on the propagation di- 
rection, but it does not experience frequency changes. The net result of this is 
that the phase of the electromagnetic wave differs from direction to direction of 
the xy plane, and this makes a GW amenable to detection by interferometric 
principles. 
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Looked at in this way, the masses represented in figure ^ only play a passive 
role in the detector, in the sense that they simply make the interference between 
the two light beams possible by providing physical support for the mirrors and 
beam splitters. In other words, the physical principles underlying the working 
of an interferometric GW detector must have to do with the interaction between 
GWs and electromagnetic waves rather than with geodesic deviations of the 
masses. 

Admittedly, this is not the most common point of view |21|. It can however 
be made precise by the following considerations, which are studied in depth in 
references 01 and f[23|. 



4.1 Test light beams in a GW- warped space-time 

According to the above, it appears that we must address the question of how 
GW-induced fluctuations in the geometry of a background space-time affect the 
properties (amplitude and phase) of a plane electromagnetic wave - a light beam - 
which travels through such space-time. It will be sufficient to consider that the 
light beam is a test beam, i.e., its back action on the surrounding geometry is 
negligible. 

Let then A^(x, t) be the vector potential which describes an electromagnetic 
wave travelling in vacuum; A p thus satisfies Maxwell's equations 

CU M = , (17) 
where □ stands for the generalised d'Alembert operator: 

□ EE gP° V p V CT . (18) 

We need only retain first order terms in h pl , in the covariant derivatives, so 
that equation ( |l7| ) reads 



□ A„ ~ rj pa (d P d a A^ - 2^ d p A v - r v pa d v A^ - A v d p r; a ) = , (19) 

To this equation, gauge conditions must be added. We shall conventionally adopt 
the usual Lorentz conditions, V^A^ = 0, which, to lowest order in the gravita- 
tional perturbations, read 

V^A" ~ -d t A t + d x A x + d y A y + d z A z - h» v d„A v = . (20) 

In addition to the weakness of the GW perturbation, it is also the case in 
actual practice that: 

• The GW typical frequencies, to, are much smaller than the frequency of the 
light, Q: lu < Q. 

• The wavelength of the GW is much larger than the cross sectional dimensions 
of the light beam. 
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Wave front distortions and beam curvature are effects which can also be safely 
neglected in first order calculations ^3|. Finally, I shall make the simplifying 
assumption of perpendicular incidence, i.e., the incoming GW propagates in a 
direction orthogonal to the interferometer's armsQ. 

We shall thus consider one of the interferometer arms aligned with the x- 
direction, and the other with the y-direction, while the incoming GW will be 
assumed to approach the detector down the z — axis. The interaction GW- light 
beam will thus occur in the z = plane, hence the GW perturbations can be 
suitably described by a function of time alone, i.e., 

h^{x,t) — ► h^ujt) , (21) 

where u> is the frequency of the GW, which we can also safely assume to be 
plane-fronted, since its source will in all cases of interest be far removed from 
the observatory. In addition, for a beam running along the x — axis, the elec- 
tromagnetic vector potential will only depend on the space variable x and on 
time, i.e., 

Ap(x,t) — ► Ap(x,t) . (22) 
The following ansatz suggests itself for a solution to equations (pj|): 



A t (t,x)=e {ujt) , (23a) 

A x (t,x) = ei (art) , (23b) 

Ay fax) =a 2 cxp[in{t-x)]e irt > 2 ^ , (23c) 

A z (t,x) =a 3 exp[itf (t-z)]^ 3 ^ , (23d) 

where eq{oj£), ei(u>t), 02 (wt) and 03(wi) are small quantities of order h, a 2 and 
<i3 are constants, and £1 is the frequency of the light. Clearly, these expressions 
reproduce the plane wave solutions to vacuum Maxwell's equations in the limit 
of flat space-time, i.e., when = 0. 

Let us now take an incoming GW of the form 

h {+>x} (t) = H {+iX} e iut , (24) 



and substitute it into (19) and (|20|), with the ansatz ( |23a| -d), neglecting higher 
order terms in the ratio ui/f2. Then p2[ , both (f>2 and 03 are seen to satisfy the 
approximate differential equation 

4 > {t) + —^{t) + 1 ^h + {t)=Q, (25) 

The solution to this equation which is independent of the initial conditions 
is, to the stated level of accuracy, 

4 This condition can be easily relaxed, but it complicates the equations to an extent 
which is inconvenient for the purposes of the present review. Details are fully given 
in reference teal. 
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n - 



2lo 



H + sin cot . 



(26) 



As shown in reference p2| , we need not worry about either £o( w O or ^i(^) 
at this stage because the longitudinal component of the electric field (i.e., E x ) 
is an order of approximation smaller than the transverse components, which are 
given by 



hence 



E y — d t A y 



E,. = 



(27) 



Ey{x,t) ~ i/2a2 exp 
E z (x,t) ~ \fi a% exp 



. . J? - 
iJ?(i — x) + i — iij- svciuot 
2to 

i£2(t — x) + i — iS+ sinwt 
2uj 



(28a) 
(28b) 



These expressions beautifully show how the incoming GW causes a phase 
shift in an electromagnetic beam of light. Note that this phase shift is a periodic 
function of time, with the frequency of the GW. If we consider a real interfer- 
ometer, such as very schematically shown in figure |^, and call r the round trip 
time for the light to go from Mo to M2 and back, then the accumulated phase 
shift is, according to these formulas, 



n - 

S x (j> = 2 x — # 4 
2lo 



sin ■ 



(29) 



since there is an obvious symmetry between light rays travelling to the right and 
to the left for a GW arriving perpendicu larly to them0. 

The arguments leading to equations (28a-b) can be very easily reproduced, 
mutatis mutandi, to obtain the phase shift experienced by a light ray travelling 
in the y, rather than the x direction - everything in fact amounts to a simple 
interchange {x < — ► y} in the equations, which includes {h + < — > —h+} as this 
is equivalent to {h xx < — ► h yy }, see (|l2|). The result is 



Sy4> 



-2 x 



ft ~ LOT 

H+ sin ■ 



(30) 



2co 2 

In the actual interferometer, provided it has equal arm lengths, the two laser 
rays recombine in the beam splitter with a net phase difference 



= 5 x (j) - S y 4> = - — H+ sin ^ 

LO 2 



and this produces an interference signal, which is in principle measurable 
the instrumentation is sufficiently sensitive. 



(31) 
- if 



5 The reader is warned that this symmetry does not happen if the GW and the light 
beam are not perpendicular, see ]23[|. 
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The reader may wonder how is it that the detector signal only depends on 
one of the GW amplitudes, h+, but not on the other, h x . The reason is that 
we have made a very special assumption regarding the orientation of the GW's 
polarisation axes relative to the light beam propagation directions. In a realistic 
case, even if perpendicular GW incidence happens, the detector's arms will not 
be aligned with the GW's natural axes, let alone the most likely case of oblique 
incidence. An important conclusion one should draw from this section is a con- 
ceptual one, that interferometric detectors are able to measure GW amplitudes 
and polarisations as a result of the interaction between the electromagnetic field 
of light rays and the background space-time geometry they travel across. 

Beyond this, though, equation ( |3l|) has very relevant quantitative conse- 
quences, too. For example, as stressed in reference its range of validity 
is not limited to interferometer arms short compared to the GW wavelength. 
Therefore, according to the formula, a null effect (signal cancellation) happens 
if the round trip time r equals the period of the GW, 2tt/uj. Likewise, equa- 
tion ( |3l| ) also tells us that maximum detector signal occurs when r = ir/ui. All 
this happen to be true for arbitrary incidence and polarisation of the incoming 
GWs as well. For GW frequencies in the 1 kHz range, the best detector should 
thus have arm lengths in the range of 150 kilometers - and even longer for lower 
GW frequencies. No ground based GW antenna has ever been conceived of such 
dimensions yet there are intelligent ways to store the light in shorter arms for 
suitably tuned GW periods. I shall not go into details of these technical matters, 
see [|| and [^lj for thorough information. 

5 Acoustic GW detectors 

Acoustic GW detectors work based on a completely different concept - see fig- 
ure 0: the idea is to set up test masses M linked together by a spring of relaxed 




U W -I 

Fig. 4. The acoustic detector concept: a GW coming perpendicular to the spring drives 
the two masses M at its ends to oscillate at the GW's frequency. Resonant amplification 
is obtained when the latter equals the spring's characteristic frequency, ft. 
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length £q, so that GW tides drive their oscillations around the equilibrium po- 
sition, with significant mechanical amplification at the spring's characteristic 
frequency fl. The spring deformation 

?(f) = £(t) - £ (32) 
thus obeys the following equation of motior|^ 

0(f) + 2 q{t) = \t a h<t) , (33) 

where 



h(t) = [h+ (f ) cos 2^ + h x (t) sin 2<p] sin 2 9 , (34) 

as follows from ([13]) - see also reference [Q for a complete discussion of this 
case. 

This is the main idea, but in practice acoustic GW detectors are elastic solids 
rather than a single spring, i.e., they do not have a single characteristic frequency 
but a whole spectrum. The response of an elastic solid to a GW excitation is 
assessed by means of the classical theory of Elasticity, as described for example 
in . In such theory the deformations of the solid are given by the values of a 
vector field of displacements, u(x, f), which satisfies the evolution equations 

Q^-^ 2 u-(X + n)V(V-u)=f(x,t), (35) 

where g is the (undeformed) solid's density, and A and /i are its Lame coefficients, 
related to the Poisson ratio and Young modulus of the material the solid is 
made of |25|| . The function in the rhs of the equation is the density of external 
forces driving the motion of the system; in the present case, these are the tides 
generated by the sweeping GW, i.e., 



fi(x,t) = gRutj(t)xj , (36) 

where Rtuj (f ) are components of the Riemann tensor evaluated at a fixed point 
of the soli d, m ost expediently chosen at its center of mass. As already discussed 
in section 2.2, plane GWs have at most six degrees of freedom, adequately asso- 
ciated with the six electric components of the Riemann tensor, Rut jit)- The six 
components are one monopole amplitude and five quadrupole amplitudes, and 
this important structure is made clear by the following expression of the density 
of GW tidal forces: 



f(x,t) = f^(x)g^(t) + J2 f (2 ' m) (x)9 (2 ' m) (t), (37) 

m= — 2 



6 I shall not include any dissipative terms at this stage, for they do not influence the 
key points of our present discussion. 
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which is entirely equivalent to ( pq ) - see reference |2q] -, and uses the common 
(l,m) notation convention to indicate multipole terms. It is very important to 
stress at this stage that f( l > m \x) are pure form factors, simply depending on 
the fact that tides are monopole-quadrupole quantities, while all the relevant 
dynamic information carried by the GW is encoded in the time dependent co- 
efficients g( l ' m \t). According to these considerations, we see that the ultimate 
objective of an acoustic GW antenna is to produce values of g( l ' m \t) - or indeed 
to extract from the device's readout as much information as possible about those 
quantities. 

Somewhat lengthy algebra permits to write down a formal solution to equa- 
tions ( |36| ) and ( |3^ ) in terms of an orthogonal series expansion p6f : 

u{x,t) = u n(x) 

N 

where 



1=0 and 2 
tn=—l,...,l 



M 



N 



I [1,1 
9n 



(38) 



/#< m) ^i- / u* N (x). (x)d 3 x, (39a) 

M J Solid 

g^ rn \t)^ f g^ m \t')smuj N (t-t')dt' , (39b) 
Jo 

with M the whole solid's mass; omega^ is the (possibly degenerate) character- 
istic frequency of the elastic body, and u^{x) the corresponding wavefunction, 
both determined by the solution to the eigenvalue problem 

H\7 2 u N + (A + fi) V(V-wat) = -uj% gu N , (40) 

with the boundary conditions that the surface of the solid be free of any tensions 
and/or tractions - see |^6| for full details. 

Historically, the first GW antennas were Weber's elastic cylinders but 
more recently, spherical detectors have been seriously considered for the next 
generation of GW antennas, as they show a number of important advantages 
over cylinders. I shall devote the next sections to a discussion of both types of 
systems, though clear priority will be given to spheres, due to their much richer 
capabilities and theoretical interest. 

5.1 Cylinders 

First thing to study the response of an elastic solid to an incoming GW is, as we 
have just seen, to determine its characteristic oscillation modes, i.e., its frequency 
spectrum lun and associated wavef unctions, un(x). In the case of a cylinder this 
is a formidable task; although its formal solution is known |p7| , ^8| , cylindrical 
antennas happen in practice tobe narrow and long p9|30[|, and so approximate 
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Fig. 5. The first longitudinal oscillation mode of a long cylinder. Note that it has a 
node at the center and maxima at the end faces. A whole period is represented, and 
transverse deformations are also shown. 



solutions can be used instead which are much simpler to handle, and sufficiently 
good - see also . 

It appears that, in the long rod approximation, the most efficiently coupled 
modes are the longitudinal ones, and these have typical sinusoidal profiles, of 
the type 

... . /nirz\ . ( nirv s t\ 

oz(z,t) oc sin y — — j sm I — — — I , n= 1,2,3,... (41) 

for a rod of length L whose end faces are at z — ±L/2, and in which the speed 
of sound is v s . Figure || graphically shows the longitudinal deformations of the 
cylinder which correspond to (fi"i"|), including transverse distortions which, though 
not reflected in the simplified equation (^Tj) , do happen in practice as a result of 
the Poisson ratio being different from zero j25|. An important detail to keep in 
mind is that odd n modes have maximum displacements at the end faces, while 
even n modes have nodes there. In fact, the latter do not couple to GWs It 
is also interesting to stress that the center of the cylinder is always a node - this 
is relevant e.g. to suspension design issues |$0| . 

A very useful concept to characterise the sensitivity of an acoustic antenna 
is its cross section for the absorption of GW energy. If an incoming GW flux 
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density of <P(l>) watts per square metre and Hertz sweeps the cylinder and sets it 
to oscillate with energy E(lo) joules then the cross section is defined by the ratio 

'-<»> = fj$ ■ m 

which is thus measured in m 2 Hz. Simple calculations show that, for optimal 
antenna orientation (perpendicular to the GW incidence direction) this quantity 



is given by 24 



8 GMv 2 s 

CTabs(^n) = o 3 , 71=1,3,5,..., (43) 

irn z c a 

where M is the cylinder's total mass. 

It is interesting to get a flavour of the order of magnitude this quantity has: 
consider a cylinder of A15056 (an aluminum alloy, for which v s ~ 5400 m/s), 
3 metres long and 60 centimetres across, which has a mass of 2.3 tonsQ; the 
above formula tells us that, for the first mode (n — 1), 



o-abs(wi) = 4.3 x 10~ 21 cm 2 Hz . (44) 

This is a very small number indeed, and gives an idea of how weak the 
coupling between GWs and matter is. 

The weakness of the coupling gives an indication of how difficult it is to detect 
GWs. By the same token, though, GWs are very weakly damped as they travel 
through matter, which means they can produce information about otherwise 
invisible regions, such as the interior of a supernova, or even the big bang. 

Equation ( (43] ) is only valid for perpendicular GW incidence. If incidence is 
instead oblique then a significant damping factor of sin 4 comes in, where 9 
is the angle between the GW direction and the cylinder's axis [^4)0- This is a 
severe penalty, and it also happens in interferometric detectors not optimally 
oriented f^|il|. 



5.2 Solid spheres 

The first initiatives to construct and operate GW detectors are due to J. Weber, 
who decided to use elastic cylinders. This philosophy and practice has survived 
himQ and still today (February 2002) all GW detectors in continuous data tak- 
ing regimes are actually Weber bars, though with significant sensitivity improve- 
ments [gU derived, amongst other, from ultracryogenic and SQUID techniques. 

About ten years after Weber began his research, R. Forward published an 
article |3^] where he pondered in a semi-quantitative way the potential virtues 
of a spherical, rather than cylindrical GW detector. Ashby and Dreitlein [Q 

7 These figures correspond to a real antenna, see 

8 Such factor can incidentally be inferred easily from (|l3|), if one notices that the 
energy of oscillation appearing in the numerator of ( ff3| ) is proportional to l 2 . 

9 Professor Joseph Weber died on 30th September 2000 at the age of 81. 
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estimated how the whole Earth, as an auto- gravitating system, responds to GWs 
bathing it, and later Wagoner Q developed a theoretical model to study the 
response of an elastic sphere to GW excitations. 

Interest in this new theoretical concept then waned to eventually re-emerge 
in the 1990s. The ALLEGRO detector group at Louisiana constructed a room 
temperature prototype antenna |36|,[37j , which produced sound experimental ev- 
idence that it is actually possible to have a working system capable of making 
multimode measurements - I'll come to this in detail shortly -, thus proving 
that a full-fledged spherical GW detector is within reach of current technologi- 
cal state of the art, as developed for Weber bars. It was apparently the fears to 
find unsurmountable difficulties in this problem which deterred further research 
on spherical GW antennas for years f38fl . 

In this section I will give the main principles and results of the theory of 
the spherical GW detector, based on a formalism which has already been partly 
used in section ||, and for whose complete detail the reader is referred to [g6| . 

As we have seen, first thing we need is the eigenmodes and frequency spectrum 
of the spherical solid. This is a classical problem, long known in the literature, 
the solution to which I will briefly review here, with some added emphasis on 
the issues of our present concern. 

The oscillation eigenmodes of a solid elastic body fall into two families: 
spheroidal and torsional modes [Efl. Of these, only the former couple to GWs, 
while torsional modes do not couple at all (39) . Spheroidal wavefunctions have 
the analytic form 

) = A n i (r) Y lm (n) n - B ni (r)inx LY lm (n) , (45) 

where Yi m are spherical harmonics |ff"o| , n = x/R is the outward pointing nor- 
mal, L is the 'angular momentum operator', L = —ix x V, A n i(r) and B n i{r) 
are somewhat complicated combinations of spherical Bessel functions |2(|, and 
{nlm} are 'quantum numbers' which label the modes. The frequency spectrum 
appears to be composed of ascending series of multipole harmonics, uj n i, i.e., for 
each multipole value I there are an infinite number of frequency harmonics, or- 
dered by increasing values of n. For example, there are monopole frequency har- 
monics wio, u>20, etc.; then dipole frequencies wn, W2i,- ■ ■ , then quadrupole 
harmonics W12, W22, and so on. Each of these frequencies is (21 + l)-fold degen- 
erate, and this is a fundamental fact which makes of the sphere a theoretically 
ideal GW detector, as we shall shortly see. 

If the above expressions are substituted into ( |39a| -b), then into (|38|), one 
easily obtains the sphere's response function as 



n=l 



(x,t) = V — u n00 (x) g { °f ] (t) 



~ 0J n2 
n—l 



U n2m(x) g^ (t) 



(46) 
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where 



~Jd Jo A »o(r)pr 3 dr , (47a) 



1 f R 

a„ 2 = - J [An2 (r) + 3 B„ 2 (r)] p r 3 dr , 



(47b) 



^(t) = / g^ m \t') smco nl (t - t') At' . (48) 



and 



o 

The series expansion (^6|) transparently shows that only monopole and qua- 
drupole spherical modes can possibly be excited by an incoming GW. The 
monopole will of course not be excited at all if General Relativity is the true 
theory of gravitation. A spherical solid is thus seen to be the best possible shape 
for a GW detector. This is because of the optimality of the overlap coefficients 
Onl between the universal form factors f^ l,m \x) and the sphere's eigenmodes 
u n im(x), which comes about due to the clean multipole structure of the latter: 
only the I = and I = 2 spheroidal modes couple to GWs, hence all the GW energy 
is deposited into them, and only them. Any other solid's shape, e.g. a cylinder, 
has eigenmodes most of which have some amount of monopole/quadrupole pro- 



jections in the form of the coefficients (39a), and this means the incoming GW 



energy is distributed amongst many modes, thus making detection less efficient. 
We shall assess quantitatively the efficiency of the spherical detector in terms of 
cross section values below. 

But, as just stated, quadrupole modes are degenerate. More specifically, they 
are 5-fold degenerate, each degenerate wavefunction corresponding to one of the 
five integer values m can take between —2 and +2. Monopole modes are instead 
non-degenerate. Figure ^ shows the shapes of all these modes - see the 
caption to the figure for further details. 

Degeneracy is a key concept for the multimode capabilities of the spherical 



detector. For, as explicitly shown by equation (46), monopole and quadrupole 
detector modes are driven by one and five GW amplitudes, respectively, i.e., 
g(°>°\t) and g^ m \t). Therefore, if one could measure the amplitudes of these 
modes, i.e., the amplitudes of the deformations displayed in figure ^, then a 
complete deconvolution of the GW signal would be accomplished. This is a unique 
feature of the spherical antenna, which is not shared by any other GW detector: 
it enables the determination of all the GW amplitudes, not just a combination 



of them, no matter where the signal comes from. In section 5.3 below I shall give 
a more detailed review of how the multimode capability can be implemented 
in practice. 



Cross sections The general definition (42) applies in this case, too. Since 
cross section is a frequency dependent concept, and since quadrupole modes 
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1 = 2 
m = +1 





:2 
: +2 




Fig. 6. The five quadrupole spheroidal oscillation modes of a solid elastic sphere, and 
the monopole mode. Note the latter (top left) is a spherically symmetric 'breathing' 
mode, while the rest have non-symmetric shapes. Because the eigenmodes ( ff5| ) are ac- 
tually complex (the spherical harmonics for m^O are complex) suitable combinations 
of them and their complex conjugates have been used to make plotting possible. These 
shapes are shared by all corresponding harmonics, i.e., quadrupole profiles for example 
are those shown no matter the harmonic number n of their frequency u n 2'. simply, they 
oscillate faster for higher n, but always keep the represented profile. 



are degenerate, it is clear that energies deposited in each of the five degenerate 
modes of a given frequency harmonic must be added up to obtain E(uj n 2) for 
that mode. Such energy must be calculated by means of volume integrals - to 
add up the energies of all differentials of mass throughout the solid -, the details 
of which I omit here. The final result turns out to be a remarkable one |Bg| : cross 
sections factorise in the form 



GM v 2 

cr a bs(uJni) = ICiC&) — - (k n ia n i) 2 (I = or 2) 



(49) 
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Table 1. Compared characteristics and cross sections for a cylindrical and a spherical 
GW detector of like fundamental frequencies. Note that the cylinder is assumed to be 
optimally oriented, i.e., with its axis perpendicular to the GW incidence direction. 



Cylinder 



Sphere 



v x = 910 Hz 

L — 3.0 metres 
D = 0.6 metres 



wa = 910 Hz 
u 22 = 1747 Hz 



2R = 3.1 metres 



M c = 2.3 tons 



M s 



42 tons 



oi = 4.3xl0~ 21 cm 2 Hz 



(Optimum orientation) 



cri2 = 9.2xl0~ iu cm 2 Hz 
<722 = 3.5 xlO" 20 cm 2 Hz 

(Omnidirectional) 



where GMv^/c 3 is a characteristic of the sphere's materialQ and (k n ia n i) a 
dimensionless quantity associated with the {nl}-th frequency harmonic; finally, 
and this is the stronger theoretical point of this expression, /Q(K) is a coefficient 
which is characteristic of the underlying theory of GWs, symbolically indicated 
with H. For example, if the latter is General Relativity (GR) then 

Ko(N) = 0, 

* = GR^{ 16?r2 (50) 



/C 2 (K) 
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while if it is e.g. Brans-Dicke |0 then these expressions get slightly more com- 
plicated Hill, etc. 

Sticking to GR, a few illustrative figures are in order. They are shown in 
table [lj where a material of aluminum A15056 alloy has been chosen. It appears 
that a sphere having the same fundamental frequency (^12) as a cylinder (u\) is 
about 20 times more massive, and this results in a significant improvement in 
cross section, since it is proportional to the detector mass. A spherical detector 
is therefore almost one order of magnitude more sensitive than a cylinder in the 
same frequency band - obviously apart from the important fact that the sphere 
has isotropic sensitivity. 

But there is more to this. Table [l] also refers to the sphere's cross section 
in its second higher quadrupole harmonic frequency, V22 ~ almost twice the 

10 v t is the so called 'transverse speed of sound', and is related to the true speed of 
sound, v s , by the formula 

vt = (2 + 2ap)- 1/2 Vs , 



with op the material's Poisson ratio. 
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Fig. 7. Cross section per unit mass of a cylinder and a sphere of like fundamental 
frequencies, in units of Gv'j/c 3 . Note the (slightly) better figure (17%) for the lowest 
mode in the sphere, as well as the appreciable value in the second harmonic of the 
latter, in sharp contrast with the null coupling of this mode in the cylinder. Third 
harmonics show a considerable reduction in sensitivity. 



value of the first, v\2. It is very interesting that cross section at this second 
frequency is only 2.61 times smaller than that at the first [Q while, as stressed 
in section |2.1^ it is zero for the second mode of the cylinder. Figure |?] shows 
a plot of the cross sections per unit mass of a cylinder and a sphere of like 
fundamental frequencies. It graphically displays the numbers given in the table, 
but also shows that, even per unit mass, the sphere is a better detector than a 
cylinder - its cross section 'curve' stays above the cylinder's. In particular, the 
first quadrupole resonance turns out to have a cross section which is 1.17 times 
that of the cylinder (per unit mass, let me stress again) Q, i.e., ~ 17% better. 
This constitutes the quantitative assessment of the discussion in the paragraph 
immediately following equation (fl8|). 



5.3 The motion sensing problem 

In order to determine the actual GW induced motions of an elastic solid a motion 
sensing system must be set up. In the case of currently operating cylinders this 
is done by what is technically known as resonant transducer JIB] ]. The idea of 
such device is to couple the large oscillating cylinder mass (a few tons) to a small 
resonator (less than 1 kg) whose characteristic frequency is accurately tuned to 
that of the cylinder. The joint dynamics of the resulting system {cylinder + 
resonator} is a two-mode beat of nearby frequencies given by 

w± (l±^v 1/2 ) , (51) 
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Fig. 8. Schematic diagramme of the coupling between a solid sphere and a resonator, 
modeled as a small mass linked to a spring attached to the sphere's surface. The dashed- 
dotted arc line on the left indicates the position of the undeformed sphere's surface, 
and the solid arc its actual position. 



where ujq is the frequency of either oscillator when uncoupled to the other, and 
r\ = M r esonator /-^cylinder- The key concept of this device is the resonant energy 
transfer, which flows back and forth between cylinder and sensor with the pe- 
riod of the beat, i.e., 2n(riujo)~ 1 ^ 2 . This means that, because the sensor's mass is 
very small compared to that of the cylinder, the amplitude of its oscillations is 
enhanced by a factor of ?y -1 / 2 relative to those of the cylinder, whence a mechan- 
ical amplification factor is obtained before the sensor oscillations are converted 
to electrical signals, and further processed - see a more detailed account of these 
principles in Q . 

The same principles can certainly be applied to make resonant motion sensors 
in a spherical antenna. In this case, however, a special bonus is there, associated 
to the degeneracy of the quadrupole frequencies: because all five quadrupole 
modes oscillate with the same frequency, it is possible to attach five (or more) 
identical resonators, tuned to a given quadrupole frequency, at suitable positions 
on the sphere surface, thus taking multiple samples of the sphere's motion. This 
makes possible to retrieve the oscillation amplitudes of the five degenerate modes 
- figure ||-, and thereby of the GW quadrupole amplitudes g^' m \t), since both 
are linearly related through equation (p6|). 

A single spherical antenna can thus deconvolve completely the quadrupole 
GW signal, and do so with isotropic sky coverage. These characteristics are 
unique to the spherical detector, and they make it a theoretically superior system 
compared to either interferometers or Weber bars. In addition, a sphere can 
naturally measure the amplitude of the non-degenerate monopole mode, as it is 
conceptually simple to sense the amplitude of an isotropically breathing pattern. 

The conceptual idea of a resonant sensor is shown in figure ||, and the equa- 
tions of motion for such a system are jlTj 
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/kV 2 m + (A + fi) V(V-u) + 

j 

^M b Q 2 b [zb(t) - u b (t)] 5 {3) {x-x b )n b + f GW {x,t) , (52a) 

6=1 

z a (t) = [z a (t)-u a (t)]+£ w (t) , a=l,...,J, (52b) 

where M a and J? a are the mass and characteristic frequency of the a-th resonator, 
/gwOM) is the GW tide on the sphere - see @ -, and ^ w (t) is the GW 
induced tidal acceleration on the resonator itself, relative to the centre of the 
sphere; 5^ is the three dimensional Dirac density, i.e., point-like connections 
between sphere and sensors are assumed. The mathematical detail of the analysis 
of these equations is somewhat sophisticated. The interested reader will find 
complete information in j|7| ; the rest of this section will be devoted to a brief 
discussion of the main conclusions of that analysis. 

First thing to stress is that equations ( |52a| -b) cannot be solved analytically, 
they must instead be solved by a perturbative procedure. The small perturbation 
parameters are the ratios 

Va ^ ~W ' Va < 1 ' a=l,...,J, (53) 

where M is the sphere's total mass. Actually, the analysis assumes that the 
resonators are all identical, any deviations from this being eventually assessed by 
suitable methods The fundamental result links the spring deformations 

to the the GW amplitudes g^ l,m \t) by the following formula, expressed in terms 
of Laplace transforms - noted with a caret ("): 



d u 



q a (s) = V- 1/2 Y / A a m) ( s ' n )9 {l ' m) (s) , a=l,...,J, (54) 

where it is assumed that the resonators' frequency fi is tuned to either a mono- 
pole or a quadrupole harmonic of the sphere. The transfer functions Aa (s; fi) 
naturally depend on whether a monopole or a quadrupole mode is selected for 
resonator tuning; I will quote here only the quadrupole case, as it is the most 
interesting one - see again full information in |p7j : 



b=l {Ca^O 



s z + u. 



<:+) 



[S +U) r 



(c) (e)* 
Va V b 



Y 2m (nb) 5n , 



(55) 



The detection of Gravitational Waves 



25 



where Vu is the c-th normalised eigenvector of the matrix P2{n a -n\ ) ), associated 
to its non-zero eigenvalue Q , P2 is a Legendre polynomial, and n a is the position 
of the a-th resonator on the sphere's surface. Finally, 



' J a± 



I = \A n2 (R)\ CaV 1/2 



a = 1, . 



(56) 



Equations (|5^)-(|56|) are the key to the GW signal deconvolution problem: 
they show that beats occur around the tuned frequency (o>„2 m this case), and 
that the resonators oscillate with amplitudes enhanced by a factor rf 1 ! 2 , as 
indeed expected. Note that these beats have frequencies which depend on the 
resonators' positions n a , as shown by the presence of the £ a coefficient in ([36]). 

The deconvolution problem consists in inferring the GW amplitudes g^' m )(s) 
from the telescope's readouts q a {s). Thus at least 5 sensors must be attached 
to the sphere's surface if e.g. the 5 quadrupole amplitudes are looked for: once 
the corresponding five q a (s) are measured the system ( jsif ) is solved for g^ 2 ' m \s), 
that's it. 

Crucial at this point is to decide where to implant the resonators, as such 
decision bears fundamentally on the simplicity, or even the possibility of solving 
the posed problem. There are two major proposals in the literature for this, 
and they are displayed in figure ^. They permit the definition of so called mode 
channels, which are linear combinations of the system readouts q a {s) which are 
directly proportional to the GW amplitudes g( 2 ' m \s). They happen to be of the 
form Ml 



TIGA 



PHC 




Fig. 9. The TIGA and PHC resonator distributions. In the former, six sensors are 
attached to the pentagonal faces of a truncated icosahedron, while in the latter there 
are two sets of five quadrupole sensors, respectively tuned to the first quadrupole 
frequency (squares) and to the second (triangles); there is an 11-th sensor (star) which 
is tuned to a monopole frequency. The relevant common characteristic of these layouts 
is that they enable the definition of mode channels - see text for details. 
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5 or 6 



m 



— — 2 2 

— ; ■ • ■ i ) 



(57) 



a=l 



both for TIGA and PHC. The actual result of these linear combinations is the 
following: 



i.e., they are convolution products of the signal and the system beats. 

This formula appears to be very powerful, as it shows that suitable sensor 
systems enable a single spherical detector to fully deconvolve all the GW ampli- 
tudes. One should however be careful about this conclusion, for the formula also 
indicates that the relevant information can only be obtained at the resonance 
frequencies oj m ± - in an ideal, non-dissipative system. In a real system, resonant 
linewidths are never infinitely sharp, they have instead a certain breadth; this 
actually makes possible the observation across wider bandwidths, provided the 
amplifier noise can be kept sufficiently small - I shall briefly come to this below 
in section pi 



5.4 Hollow and dual spheres 

The real merit of the just described spherical GW detector comes from its sym- 
metry. A hollow sphere does of course share the symmetry properties of a solid 
sphere, and one might therefore expect it to be an interesting alternative, too. 
A detailed study of the performance of a hollow spherical GW detector can be 
found in f49[] . The added bonus of a hollow sphere is that there is one more 
structural parameter one can adjust to enhance this or that property, and this 
is of course the thickness of the spherical shell. It appears that cross sections at 
e.g. different quadrupole harmonics have characteristic behaviours when plotted 
as functions of thickness; also, one can decide to attach the motion sensors to 
the inner or to the outer side of the shell, as their GW induced oscillations have 
different amplitudes. 

To actually build and suspend a hollow sphere in the laboratory may be a 
difficult task from a technological point of view. Recently, though, a new concept 
spherical detector appeared in the literature |50|j : this is called dual sphere, and 
consists in a solid sphere nested inside a hollow one, concentric with it and with a 
narrow gap between them - figure [10[ An incoming GW drives both spheres into 
oscillation. Clearly, the inner solid sphere has a higher first (quadrupole) mode 
frequency than the outer hollow one, therefore an incoming GW with a frequency 




(58) 
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Fig. 10. Conceptual scheme of the dual sphere: the GW signal drives the facing sur- 
faces between the nested spheres into oscillatory motions; when the signal frequency 
falls between the resonances of the spheres, oscillations happen with phase opposition, 
thereby enhancing the device's response. Motion sensing can be accomplished e.g. by 
non-resonant optical transducers. 



between those two will drive the oscillations of the spheres with opposite phasej n |, 
thus enhancing the signal by a rough factor of 2. 

The motion sensing in either hollow or dual spheres is conceptually analogous 
to that in the solid sphere, with the added flexibility that in the hollow piece 
one can choose to sense the displacements of either side of the shell. Actually, 
though, it appears that non-resonant detectors appear constitute a better choice 
in dual spheres, for this enables a significant bandwidth enlargement (5(J. 

Let me briefly discuss now, for completeness, a few essentials of GW detector 
sensitivity. 



6 GW detector sensitivities 



So far I have only discussed the theoretical basis of the workings of GW detectors, 
whether interferometric or acoustic, yet have made no mention of the practical 
difficulties of getting them actually working. . . 

The extreme weakness of any expected GW signals arriving in the Earth J52[ 
is in fact a source of such truly difficult problems that it has prevented GW 
detectors from sighting a real signal in the last 40 years, since the times of 
J. Weber. Local detector noise has to date overwhelmed any signals possibly 
hitting the antennas, and therefore the technological challenge has been for years, 

11 As shown in textbooks on Mechanics, see e.g. fl5l[, there is phase change in an 
oscillator's response to a periodic excitation as the frequency of the latter shifts 
from below the oscillator's natural frequency {Q, say) to above it; the transition 
region has a width of order fl/Q centred at Q, where Q is the oscillator's mechanical 
quality factor, about 10 7 or more in GW detectors. 
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and still is today, to reduce that noise to the level where it can be filtered out 
with a meaningful signal to noise ratio |^3| . 

During the last decade or so, a number of people in different countries world- 
wide have managed to get important GW detection research projects funded 
which constitute a major step forward in detector technology. Their goal is 
to improve the sensitivity to the point where a significant event rate becomes 
available to the GW astronomer. We thus find such laboratories as VIRGO (a 
French- Italian collaboration), LIGO (USA alone), or LISA (a space mission, 
jointly funded by NASA and ESA, the European Space Agency). In addition to 
those, somewhat smaller experiments are GEO-600 (a German-British venture) 
and TAMA (the Japanese project, currently making strong progress in both and 
stability). 

So much for interferometric antennas. But endeavours have not declined in 
the acoustic detector arena, either. In fact, the five acoustic detectors of the 
Weber type (NAUTILUS, EXPLORER, AURIGA, ALLEGRO and NIOBE) 
constitute the only working systems in the world today. Unfortunately, though, 
they are only sensitive to catastrophic events happening in our galaxy, with a far 
too low occurrence rate. These detector systems are periodically upgraded, and 
their sensitivity has gone up a few orders of magnitude since their origins. As 
already stressed in section 0, new generation spherical GW detectors are being 
programmed in Brazil, Holland and Italy, and these should suddenly improve 
over bars by at least one more order of magnitude. 

Figure [ll] contains a recent plot of GW sensitivity of various earth based 
detectors, therefore in the frequency range near 1kHz. On the other hand, LISA 
will be sensitive at frequencies far away from the range plotted in figure [ll] - see 
LJSA's web site at http : / / sci . esa. int/home/lisa and figure 12 below. What 
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we see in ordinates in plots like the one in figure [ll| which are standard in GW 
science, is a spectral density - or, rather, its square root. It is to be understood 
as follows. 

Noise of whatever origin causes any detector to generate random outputs 
- stochastic time series, as this is technically known. This noise competes with 
any GW signals which may be present in the antenna readout, obscuring their 
detection. Noise of course does not come from the sky, it rather gets added to the 
GW signal at almost all the different stages of the detection process. For example, 
the GW induced oscillations of a given mass compete with thermal oscillations 
of that mass and its suspension systems; then there is noise in the conversion 
from mechanical or optical GW signals to electrical signals; then there is noise 
introduced by the electronics in the amplifiers of the latter; and so forth. . . 

In the end, the detector output is a certain physical magnitude (e.g. a volt- 
age) which consists of a number of various superimposed sources of noise plus, 
possibly, a GW signal converted to volts by the detector hardware. One infers 
the value of the actual dimensionless GW amplitude h(t), say, from the output 
voltage if one knows the precise physics of the transduction process - and this 
is obviously the case in any useful instrument. 

Conversely, it is also possible to back- convert all the sources of noise voltage 
picked up across the various detector stages to an equivalent dimensionless "GW 
noise", which gets directly added to real GW signals, and travels through an 
ideally noiseless detector. This artifact is expedient because, for a fixed antenna, 
it enables a quick assessment of the detectability of a given GW signal, normally 
calculated by methods of gravitation theory, by direct comparison with suitably 
constructed detector characteristic curves - such as those in figure 

Let us then call x(t) any one of the detector's readout channels, back-con- 
verted to a GW amplitude by the above described procedure. We split this up 
into a GW signal proper h(t) plus a noise term n(t): 

x(t) = h(t) + n(t) . (59) 

For stationary Gaussian noise, the statistical properties of n(t) are encoded 
in its spectral density function, Sh(w), which is the Fourier transform of the 
autocorrelation function 

/oo 
R(t) e iuJT dr , (60) 
-oo 

where (— } stands for ensemble average ]54| ]. It can be shown |5i| that the op- 
timum filter to extract the signal h(t) from the system readout x{t) is the so 
called matched filter, whose transfer function is the signal's Fourier transform 
h(u>) divided by the spectral density Sh(oj), and the detection threshold can be 
set from the integrated signal to noise ratio: 



(61) 
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Fig. 12. Root mean square noise spectral densities, referred to GW amplitudes, for 
some of the upcoming interferometric GW detectors, together with the spectral inten- 
sities of various signals. The latter are estimated by numerical calculation, while the 
noise is modeled on the basis of its origin and instrumental characteristics. See |5q]. 



The GW signal hit) is a dimensionless quantity, as it measures a perturba- 
tion of the Minkowski metric, see equation (|^). Therefore the spectral density 
Sh(uj) has dimensions of time, or inverse frequency, Hz -1 , accordin g to the def- 
initions (|6(]). Now, signal to noise ratio as denned by the integral ( pl| ) is made 
up of the contributions of the ratio between the signal power |/i(u)| 2 to the noise 
power Sh(u) at all frequencies; the idea of the graphical representation in fig- 
ure ^ is thus to show which is the level of noise at different frequencies by means 
of an rms quantity, such as is the square root of the noise spectral density. The 
appropriate units for this representation are accordingly Hz~ 1//2 , as indicated. 

Different interesting sources of GWs are being considered by other authors in 
this volume, so I will not go into such matters here. It is nevertheless instructive 
to present an example graph of a few signals on top of various detectors' sensitiv- 
ity curves in order to get a picture of the actual possibilities of each instrument, 
and also to grasp what are the spectral orders of magnitude of different GW sig- 
nals. One such plot is presented in figure This is, let me insist, the standard 
way to assess the detectability of the different GW sources. 



The detection of Gravitational Waves 



31 



7 Concluding remarks 

This article is a brief review on the nature of the interactions of GWs with test 
particles and test electromagnetic fields, as they specifically happen in currently 
conceived detection devices. While the fundamental principles are not new in 
themselves, their application in actual systems is still in many cases subject 
of research, as we have seen. A thorough understanding of these matters is 
absolutely essential for an adequate interpretation of the antenna readouts, the 
more so if one considers the extreme weakness of any signals reaching us from 
even the most intense sources. 

I have omitted any detailed discussion of the practical problems faced by 
real detector building. This is a major research field in itself, of an intrinsically 
sophisticated and multidisciplinary nature - involving such fundamental issues 
as quantum measurement limits and techniques, or Quantum Optics |pq| . But it 
is of course not directly related to Astrophysics or Relativity. . . I have however 
considered appropriate to summarily brief the interested reader on how detector 
sensitivities are defined, and on which are the detection thresholds in current 
state of the art GW detectors: this is a key issue for an astrophysicist/relativist 
wishing to assess the detectability of a given GW signal, whether by existing 
instruments or future planned. 

GW detection endeavours have been the subject of intensive research during 
the last 4 decades, though the ultimate goal of sighting GW events has still not 
been accomplished. While this may look like a major failure, one may not forget 
that detector sensitivities have gone up by a remarkable six orders of magnitude 
(in energy, three in GW amplitude) since the very first telescopes constructed 
by J. Weber. Also, a look at the trend in this development indicates that we are 
getting closer to the obejctive. 

Prospects look now better than ever yet the real challenge is still there. . . 
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